How to compute one-loop Feynman diagrams in lattice QCD 
with Wilson fermions 
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We describe an algebraic algorithm which allows to express every one-loop lattice integral with gluon or Wilson- 
fermion propagators in terms of a small number of basic constants which can be computed with arbitrary high 
precision. Although the presentation is restricted to four dimensions the technique can be generalized to every 
space dimension. We also give a method to express the lattice free propagator for Wilson fermions in coordinate 
space as a linear function of its values in eight points near the origin. This is an essential step in order to apply 
the recent methods of Liischer and Weisz to higher-loop integrals with fermions. 



In [y we presented a general technique which 
allows to express every one-loop bosonic integral 
at zero external momentum in terms of two un- 
known basic quantities which could be computed 
numerically with high precision. This method al- 
lows the complete evaluation of every diagram 
with gluon propagators. 

We have recently generalized [|| this technique 
to deal with integrals with both gluonic and 
fermionic propagators. For the fermions we use 
the Wilson action [||. Notice that our method 
depends only on the structure of the propagator 
and thus it can be applied in calculations with 
the standard Wilson action as well as with the 
improved clover action [Q. We show that every 
integral at zero external momentum can be ex- 
pressed in terms of a small number of basic quan- 
tities (nine for purely fermionic integrals, fifteen 
for integrals with bosonic and fermionic propaga- 
tors). The advantage of this procedure is twofold: 
first of all every Feynman diagram can be com- 
puted in a completely symbolic way making it 
easier to perform checks and verify cancellations; 
moreover the basic constants can be easily com- 
puted with high precision and thus the numer- 
ical error on the final result can be reduced at 
will. Although the presentation is restricted to 
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four dimensions the technique can be generalized 
to every space dimension. 
Define 

an Hy Hf. 

^A'^^Y Dl+\k,mf)Dl{k,mb) 

where Ui are positive integers, p and q are ar- 
bitrary integers (not necessarily positive), fc^ = 
2sin(fc,,/2) 

DB{k,mb) = k^ + ml^^kl +ml; (2) 

i 

DpiKruf) = ^sin^fc, + ^(P)2+m^. (3) 

i 

is the denominator appearing in the propagator 
for Wilson fermions^. In the following when one 

^To be precise, Dp{k, mj) is the denominator in the prop- 
agator for Wilson fermions only for ruf = 0. For mf ^ 
the correct denominator would be 

Dpik, mf) = '^sm'^ (j^P +mf^ . (4) 

i 

However in our discussion m j will only play the role of an 
infrared regulator and thus it docs not need to be the true 
fermion mass. The definition (H) is easier to handle than 
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of the arguments Ui is zero it will be omitted as 
an argument of J-". The parameter S is used in 
the intermediate steps of the calculation and will 
be set to zero at the end. 

To simplify the discussion we have only consid- 
ered the case rw = 1 but the technique can be 
applied to every value of rw- Moreover we have 
restricted our attention to the massless case, i.e. 
we have considered the integrals J-s in the limit 
rub = nif = m — > 0. 

We have developed a procedure that allows to 
compute iteratively a generic which is J-'g in 
the limit ^ — > 0, in terms of a finite number of 
them: precisely every J-{p,q;nx,ny,nz^nt) with 
p > and q < (purely fermionic integrals) 
can be expressed in terms of jr(l,0), 
H^.~2), ^(2,0), .F(2,-l), T{2,-2), T{i,-2), 
J^(3,-3) and J^(3,-4); the integral T{2,0) ap- 
pears only in infrared-divergent integrals and we 
write it as 

^(2,0)= -^(logm2+7£-Fo) +yo (5) 

where lo is a numerical constant. If g > the 
result contains three additional constants which 
we have chosen to be 

Yi = i.F(l,l; 1,1,1) 

>2 = ^.F(l,l;l,l,l,l) 

^3 = ^.F(l,2;l,l,l) 



(6) 



together with the bosonic quantities Zo, Zi and 

Fo - "fE- 

As a first step in our procedure we express each 
integral J-sip, q;nx,ny,nz,nt) in terms of inte- 
grals of the form !Fs{r,s) only. This is achieved 
by the use of the identity 



4 
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in the case one of the n^'s is 1. In the case one of 
the rii's is 2 we use instead 

4 

= 4.[DB{k,m)~DF{k,m)] 

1=1 

+ [DB{k,m)-rr?]^ (8) 



In all other cases we use an integration by parts. 
It is here that we found convenient to have a non- 
vanishing 5, otherwise integration by parts in the 
case in which there is only one fermionic propa- 
gator would produce a \ogDp. 
As a result we have that 



Ts{p,q;nx,ny,nz,nt) = 

p q+k 
r—p~k-\-l s—q—k 



(9) 



where k = [ux + Uy + Uz + rit), and arsim, S) is a 
polynomial in m^. 

At the second step we express every Tsip, q) in 
terms of a finite number of them. We start by 
inserting the trivial identity 



= 1- 



(10) 



Ds(fc,m) 

into the integral J-s{p, q; 1, 1, 1, 1) to get 

O^Ii{p,q) = 

Tsip, q; 1, 1, 1, 1) - 4Ts{p, g + 1; 2, 1, 1, 1) 
-m^Ts{p,q+l;lA,l,l) (11) 



Using the previous recursions we obtain from (|11|) 
a non trivial relation among the .^-^(r, s). Start- 
ing from the analogous relation obtained from the 
fermionic propagator we get a new set of identi- 
ties I2 {p, q) which can be used to provide addi- 
tional relations between the remaining integrals. 
We end up with the result we have quoted above. 

In order to numerically evaluate the eight ba- 
sic purely-fermionic integrals we considered the 
integrals 



Jq = -q) 6 < q < 13 



(12) 



which can be analytically expressed in terms of 
the elements in the fermionic basis. We computed 



1 ^ JB(fc,0)g 

^4 2^ DF{k,0) 



(13) 



where k runs over the points k ~ {2'k/ L)(ni + 
^, 712 -I- 5, ^3 -f i, 714 + i), < 77i < L for various 
values of L between 50 and 100 and extrapolated 
the results to infinite volume by using 



Jq,L ~ Jq [1+ j^2q+2 



(14) 
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The inversion of this set of equations provides a 
very accurate estimate of the elements in the ba- 
sis (see Table 1). Analogously to compute lo, Fi, 
Y2 and Y3, we computed numerically ^(1,1,8), 
J^(2,l,9), J^(3, 1,10) and JP(5,2,11) and then 
solved the corresponding equations. 

A second important application of our method 
is connected with the use of coordinate-space 
methods for the evaluation of higher-loop Feyn- 
man diagrams. This technique, introduced by 
Liischer and Weisz [|| , is extremely powerful and 
allows a very precise determination of two- and 
higher-loop integrals. One of the basic ingredi- 
ents of this method is the computation of the free 
propagator in coordinate space. But this can be 
achieved by the use of our relations. Consider 



G{p, q; x) 



dk 



(27r)4 DP,(k,m)D%{k,m) 

dk Dp cos kf^Xf, 
(27r)4 D^p{k,m)D%{k,m) 



(15) 
(16) 

Then express cos(fc^x^) as a polynomial in k'^. It 
follows that G is a linear combination of the J- 
and therefore they can be expressed on our basis. 

Here, like in the bosonic case, the expression of 
G in terms of the basis become numerically un- 
stable for I a; I — > 00: numerical errors in the basis 
become amplified! This problem has a standard 
way out: if the expressions are unstable going out- 
ward from the origin, they will be stable in the 
opposite direction: thus, if we want to compute 
the propagator for |a;[ < d, for some fixed d, we 
choose eight points with \x\ « d (say yi, . . . , yg) 
and then we express the propagator for < d in 
terms of G(p, g; j/i), i = 1, . . . , 8. The new expres- 
sions are numerically stable: the numerical error 
on G(j),q;yi) gets reduced when we compute the 
propagator for 0. As noticed in § the 

instability of the recursion can also be used to 
provide precise estimates for the basic constants. 
We have thus used this method to obtain an in- 
dependent numerical estimate of the eight purely 
fermionic and infrared finite constants in our ba- 
sis, considering the set of eight points X*^"^ = 
{(n, [0-3],0,0),(n+l,[0-3],0,0)}. 

Because of better convergence, by using larger 
negative values of q, one can obtain more pre- 



^(1,0) 


0.08539036359532067913516702879 
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^(1,-1) 


0.46936331002699614475347539703 


^(1,-2) 


3.39456907367713000586008689687 
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0.05188019503901136636490228763 
0.05188019503901136636490228706 


m-2) 


0.23874773756341478520233613924 
0.23874773756341478520233613767 


^(3, -2) 


0.03447644143803223145396188143 
0.03447644143803223145396188128 


^(3, -3) 


0.13202727122781293085314731096 
0.13202727122781293085314731035 


^(3, -4) 


0.75167199030295682253543148585 
0.75167199030295682253543148379 



Table 1 

Numerical values of the constants appearing in 
the fermionic integrals evaluated by the two dif- 
ferent methods described. The first line refers to 
the method which goes through the evaluation 
of the integrals Jq, the second one goes through 
the propagators in coordinate space. In all cases 
there is an agreement of at least 27 digits. 



cise estimates of G(l,(7;x) at the set of points 
Evaluating the integrals by computing dis- 
crete sums on finite lattices, as before, with sizes 
L = 50 - 100, we obtain the values of G(l, —3; x) 
at the points and from these numbers we 

obtained an independent determination of the fi- 
nite eight fermionic constants (see Table I). 
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